Strong-Coupling Expansions for Multiparticle Excitations: 
Continuum and Bound States 
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We present a new linked cluster expansion for calculating properties of multiparticle excitation 
spectra to high orders. We use it to obtain the two-particle spectra for systems of coupled spin-half 
dimers. We find that even for weakly coupled dimers the spectrum is very rich, consisting of many 
bound states. The number of bound states depends on both geometry of coupling and frustration. 
Many of the bound states can only be seen by going to sufficiently high orders in the perturbation 
theory, showing the extended character of the pair-attraction. 
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The study of bound states, life times and spectral 
weights for multiparticle excitations remains a challeng- 
ing problem in many-body physics. From an exper- 
imental point of view, a variety of probes on several 
low-dimensional magnetic and strongly correlated elec- 
tronic systems show spectral features associated with 
multiparticle continuum and bound states. These fea- 
tures in two-magnon Raman spectra, optical absorption 
spectra, photoemission and even in neutron scattering 
spectra remain poorly understood. From a theoretical 
point of view, one of the most intriguing issues is the 
role the growing number of bound states play in the 
confincmcnt-dcconhncment transition in spin-Pcierls sys- 
tems as the spectrum completely changes from soliton- 
antisoliton continuum to triplets, their bound states and 
continuum Existing computational approaches to 

these problems are not adequate. 

In this letter, we present significant advances in calcu- 
lating multi-particle spectral properties from high-order 
strong coupling expansions. At the heart of our method is 
a generalization of Gclfand's linked cluster expansion for 
single-particle excited states to multi-particle states. 
From a technical point of view, our most notable achieve- 
ment is the development of an orthogonality transforma- 
tion which leads to a linked cluster theorem for multi- 
particle states even with quantum numbers identical to 
the ground state. Our method is quite distinct from the 
flow equation method of Wegner , which has been used 
recently by Knetter et al. Q for also studying multi- 
particle spectral properties, but for a more restricted 
class of models with equispaced unperturbed eigenvalues 
and an exact direct-product ground state. 

We apply the method to systems of coupled spin-half 
dimers in various one and two-dimensional geometries. 
The calculations greatly simplify when we deal with mod- 
els where the ground state is known exactly to be a prod- 
uct of dimers, such as the Shastry-Sutherland models in 
one and two dimensions 0]. Here, we focus attention on 



the one-dimensional systems. Results are presented for 
a spin-ladder system, the alternating spin-chain model, 
and the Shastry-Sutherland model. 

Except for the spin-ladder model, the high order calcu- 
lations produce qualitatively unexpected results. There 
are multiple bound states with 5 = 0,1 and 2 even in the 
weakly-coupled dimer limit. The number of bound states 
varies with frustration and many of them only show up 
when the expansions are done to sufficiently high order. 
This presumably is due to the extended range of the 
attractive interaction needed to see the multiple bound 
states. The interval of k-values where the bound states 
exist also changes with the expansion parameter. 

We begin describing our method with a Hamiltonian 



H — Hq + XHi , 



(1) 



where the unperturbed part Hq is exactly solvable, and 
A is the perturbation parameter. The aim is to calculate 
perturbation series in A for the eigenvalues of H and other 
quantities of interest. As is well known the ground 
state energy and correlation functions have a 'cluster ad- 
dition property' and hence can be calculated by linked 
cluster expansion. 

We wish to consider the excited-state many-particle 
sectors of the Hilbert space, where a "particle" may refer 
to a lattice fermion, a spin-flip, or other excitation, de- 
pending on the model at hand. The key step is to 'block 
diagonalize' the Hamiltonian on any finite cluster to form 
an effective Hamiltonian, via an orthogonal transforma- 
tion (here we will only consider real Hamiltonians): 



H eS = T HO 



(2) 



where O = e and S is real, antisymmetric. This trans- 
formation is constructed order-by-order in perturbation 
theory, so that the ground state sits in a block by it- 
self, the 1-particle states (which form a degenerate man- 
ifold under Hq, in general) form another block, the 2- 
particle states another block, and so on. The off-diagonal 
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blocks of S are determined by the requirement that the 
off-diagonal blocks of H cS vanish; and we choose the di- 
agonal blocks of S to be zero. 

Let the matrix element of H cS between initial 1- 
particle state |i) and final 1-particle state labelled 
according to their positions on the lattice, be 

£i(i,j) = <j|ff cff |i). (3) 

This excited state energy is not extensive, and does not 
obey the 'cluster addition property'. However, as shown 
by Gelfand [Q, the "irreducible" 1-particle matrix ele- 
ment 

Ai(i,j)=£?i(i,j)-^j (4) 

does have the 'cluster addition property'. Furthermore, 
for a translationally invariant system, the 1-particle 
states are eigenstates of momentum: 

i k > = -/=E ex p( jk -j)ij> ( g ) 

j 

(where N is the number of sites in the lattice), with en- 
ergy above the ground state of 

wi(k)=5^A 1 (<5)cos(k-(5). (6) 
s 

To generalize to two-particle states, let 

£ 2 (i,j;k,lH(k,l|ff cff |i,j) (7) 

be the matrix element between initial 2-particle state |i, j) 
and final state |k, 1). To obtain a quantity obeying the 
cluster addition property, we must subtract the ground- 
state energy and 1-particle contributions, to form the ir- 
reducible 2-particle matrix element (see Fig. Q): 

A 2 (i, j; k, 1) = E 2 (i, j; k, 1) - E Q (6 iM 6j,i + tf u <5j, k ) 
-Aiftkjaj.i-Ax&ljtfj.k 
-Aia,k)5 M -AiCj,l)5i, k (8) 

This quantity is easily found to be zero for any cluster 
unless i, j, k and 1 are all included in that cluster, and 
it obeys the cluster addition property. The block diago- 
nalization ensures that two particles cannot "annihilate" 
from one cluster and "reappear" on another disconnected 
one. Thus the matrix elements of A 2 can be expanded 
in terms of connected clusters alone, which are rooted or 
connected to all four positions i, j, k, 1. 

Once the effective two-particle Hamiltonian is known, 
we still have to solve the Schrodinger equation. Consider 
for simplicity a 2-particle state of identical particles, or 
one symmetric under particle exchange. Then expand 
the 2-particle eigenstate 

M = X>|iJ>, = ( Q ) 

i>j 




FIG. 1. Decomposition of a 2-particle matrix element into 
irreducible components. The round box denotes the full ma- 
trix element, the square boxes the irreducible matrix ele- 
ments, and the single line denotes the delta function. 



the Schrodinger equation takes the form 
(E-Eo)fis - ^[A 1 (k,i)/ kj + A 1 (k,j)/ ik ] 

k 

= ^A 2 (k,l;i,j)/ kl 

k,l 

-A 1 (j,i)/ jj -A 1 (i,j)/ ii , alii,]. (10) 

The fictitious amplitudes fa are defined by these equa- 
tions H , and are introduced to simplify the Fourier trans- 
form. Defining the centre-of-mass momentum K and the 
relative momentum q, this can be turned into an integral 
equation: 

[E - E - 2 A i GO C0S ( K ' S/2) cos(q • S)]f(K, q) = 

s 

l^/(K,q')[i Yl A 2 (r, < 5 1 ,<5 2 )cos(K-r)cos(q- ( 5 1 ) 

q' r,<5i,(52 

x cos(q' • 6 2 ) - Ai(<$) cos(K • 6/2) cos(q ■ 6)], (11) 
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which can be solved by standard numerical techniques 
like discretization. The spectrum of the discretized 
Hamiltonian can be obtained by diagonalization. In this 
way the spectrum and even the two particle density of 
states can be obtained. The continuum is limited by the 
the maximum (minimum) of the energy of two single par- 
ticle excitations whose combined momentum is the cen- 
ter of mass momentum, which serves as an independent 
check. 

We have used this method to investigate the low- 
energy excitation spectrum of the 2-leg spin-| Heisenberg 
ladder 

H = J2 [JSi ■ Si+i + JS' t ■ S' i+ i + J±Si ■ S'i] , (12) 

i 

where the interactions along the ladder ( J) and along the 
rungs (Jj_) are assumed to be antiferromagnetic. 

For J/J± < oo, the ground state of this model evolves 
smoothly from a product of singlet states along the rungs 
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FIG. 2. Low-energy excitation spectrum of the Heisenberg 
spin ladder for J/J± = 1/2. Beside the two-particle con- 
tinuum (gray shaded) and the elementary triplet excitation 
(dotted line) there are three massive quasiparticles: a singlet 
bound state (solid line), a triplet bound state (dashed line) 
and a quintet antibound state (dashed-dotted line). 
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FIG. 3. The binding energy Ei for two singlet (Si and S%) 
and two triplet (Ti and T2) bound states versus momentum 
K for the Ji — J2 — S chain with 5 = 0.6 and a — 0. The inset 
enlarges the region near K — n/2. 



of the ladder and has a gapped excitation spectrum 
|p[-|lT|. The occurrence of two-particle bound states in 
this model has been shown by first-order strong-coupling 
expansions 0,[l|| as well as a leading order calculation 
using the analytic Brueckner approach [l^Jl5t| . 

Starting from the dimerized ground state we have cal- 
culated series in J/J± for A2 up to order 7 for singlet 
states, and to order 12 for triplet and quintet states Jl6| . 
Fig. H shows the generic shape of the two-particle contin- 
uum as well as the low-lying massive excitations. Beside 
the elementary triplet excitation the spectrum shows ad- 
ditional singlet (S — 0) and triplet (S = 1) excitations 
which are bound states of two elementary triplets. In 
the vicinity of K w 7r there is also an S = 2 antibound 
state above the continuum. At J/J± = 1/2, we find the 
binding energy for the singlet bound state at K = it is 
Eb/J± = 0.51, substantially larger than the value 0.35 
obtained in jlijj ]. 

Further, we have studied the occurrence of bound 
states in the alternating Heisenberg chain 

H = J £ K 1 + (-!)^) S * ' S *+i + ■ S *+2] > (13) 

i 

where the are again spin-^ operators at site i, a pa- 
rameterizes a next-nearest neighbor coupling and 6 is the 
alternating dimcrization. 

The two-particle excitations have been discussed in a 
leading order Brueckner ansatz calculation Jl7],[l8| , a sec- 
ond order series expansion and an RPA study p0[ . 
With our new technique, we perform high-order series 
expansions in powers of A = (1 — $)/(! + $). Here we will 
only concentrate on the expansions for the following two 
cases: 

(1) a =0, that is, without the second neighbor inter- 



action. The series for A 2 has been computed up to order 
6 for singlet bound states, and to order 11 for triplet 
and quintet states. Here we find two singlet, S\ and S2, 
and two triplet, T\ and T2, bound states below the two- 
particle continuum. The binding energy of these bound 
states versus momentum K is given in Fig. |^ for a rather 
large dimcrization 5 = 0.6. The singlet Si exists for the 
whole range of momenta, while the triplet Ti exists only 
in a limited range of momenta and the singlet S2 and 
triplet Ti bound states occur for momenta in the vicin- 
ity of the band maximum at K = tt/2. The existence of 
the second pair of bound states has not been reported by 
previous calculations, most likely due to a limited preci- 
sion or a general incapability to deal with multiple bound 
states. The binding energy at K = ir/2 versus dimeriza- 
tion S is plotted in Fig. [l| In the limit A — > 0, the binding 
energies for Si and Ti are proportional to A, as expected, 
since the formation of these bound states is due to the at- 
traction of two triplets on neighboring sites. For S2 and 
T2, we find their binding energies are proportional to A 2 . 
This means that there is a strong enough effective attrac- 
tion between two triplets separated by a singlet dimer to 
form those bound states. 

(2) a = (1 — S)/2, that is, the expansion is along the 
disorder line where the ground states are known exactly. 
The series for A 2 has been computed up to order A 19 for 
2-particle singlet, triplet and quintet states. The two- 
particle excitation spectrum for S = 0.4 is shown in Fig. 
||. Here we find that there are three singlet and three 
triplet bound states below the two-particle continuum, 
and two quintet antibound states above the continuum. 
The energy gap at K = tt/2 for one of the singlet bound 
states, Si, is 1 + 36 exactly. In the limit A — > 0, the 
binding energies for S n and T n (n = 1,2) are of order 
A", just as for a = 0, but that for S3 and T 3 are at most 
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FIG. 4. The scaled binding energy fEt at K = n/2 versus 
dimerization S for two singlet (Si and 5*2 ) and two triplet (Ti 
and T2) states of the Ji — J2 — 5 chain with a — 0, where / 
is a scale factor. 



of order A 4 . This calculation demonstrates the power of 
the method in bringing out the complex character of the 
pair-attraction in the frustrated model. 
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FIG. 5. The excitation spectrum of the Ji — J2 — 5 chain for 
8 = 0.4 and a — (1 — <5)/2. Beside the two-particle continuum 
(gray shaded), there are three singlet bound states (Si, S2 
and S3), three triplet bound states (Ti, T2 and T3) and two 
quintet antibound states (Qi and Q2)- The inset enlarges 
the region near K — tt/2 so we can see S3 and T3 below the 
continuum. 



In conclusion, we have demonstrated a new power- 
ful approach to calculate high-order series expansion for 
quantum Hamiltonian lattice models. The application to 
the Heisenberg spin ladder and the alternating Heisen- 
berg chain has resulted in precise calculations of the low- 
lying excitation spectra of these models. For the alternat- 
ing Heisenberg chain it was shown that there are multiple 
massive singlet and triplet excitations below the contin- 
uum, which depend on frustration. 

This work was initiated at the Quantum Magnetism 
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